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Abstract 



The pair distribution function and the static structure factor are computed 

for composite fermions. Clear and robust evidence for a Ikp structure is seen 

in a range of filhng factors in the vicinity of the half-fiUed Landau level. 

Surprisingly, it is found that filled Landau levels of composite fermions, i.e. 

incompressible FQHE states, bear a stronger resemblance to a Fermi sea than 

do filled Landau levels of electrons. 
71.10.Pm,73.40.Hm 
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I. INTRODUCTION 



A central feature of the composite fermion (CF) theory |jT|-0 of the quantum Hall effect 
is that the interaction energy of electrons is transmuted into the kinetic energy of composite 
fermions. In other words, while the original problem had no kinetic energy (more precisely, a 
constant kinetic energy) but only the interaction energy, the kinetic energy is the important 
energy of composite fermions, the interaction between which is weak and can be neglected 
altogether to a good first approximation. The kinetic energy of composite fermions manifests 
itself in the form of Landau levels (LLs) of composite fermions, resulting in the fractional 
quantum Hall effect (FQHE), and, in an extreme limit, as the Fermi sea of composite 
fermions at the half filled Landau level. 

In this work, we will compute the pair distribution function and the static structure 
factor for filled Landau levels of composite fermions, corresponding to FQHE states at 
V = n/(2n +1), as well as for composite fermions at zero effective magnetic field, using 
the recently developed lowest-LL representation of the CF wave functions that enables a 
treatment of large CF systems The qualitative behavior of the results provides a striking 
evidence for Fermi-sea-like correlations in a broad range of filling factors near the half-filled 
Landau level. 

The results obtained below are valid to the extent the CF wave functions are. The CF 
wave functions have been compared extensively to the exact eigenstates, available numeri- 
cally for systems containing ~ 10 — 12 particles, and found to have close to 100% overlap 
with them It is quite plausible that the excellent agreement will persist to larger 

systems and also for other incompressible FQHE states, suggesting that the CF wave func- 
tions should be reliable up to 9/19, which is the fraction closest to 1/2 where experimental 
indications for FQHE have been observed IQ. There is sound experimental evidence that 
the CF description is valid for the compressible state at z/ = 1/2 as well [[HQ. Here, the CF 
wave functions are expected to capture the short- distance behavior of the composite-Fermi 
sea state but may not to be very useful for subtleties associated with the low-energy long- 
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distance physics; unlike the FQHE states, the composite Fermi-sea is hkely to be susceptible 
to perturbations due to the absence of a gap. Fortunately, however, we will see that it is 
not necessary to go all the way to 1/2 to see Fermi-sea-like correlations; they are manifest 
even in incompressible FQHE states away from 1/2, indicating that a composite- Fermi-sea 
description of the half-filled LL state is valid, at least as a first order description. The issues 
concerning subtle and asymptotic deviations from a perfect Fermi liquid, which can surely 
be expected given that there was no reason to expect any kind of Fermi sea in the first place, 
are beyond the scope of the present study. 

Another widely used approach for an investigation of the CF state is based on a Chern- 
Simons field theory p,|TO|JTl| , useful for a description of the low-energy long-distance physics 
of the compressible composite- Fermi sea at u = 1/2, which has been found to be reasonably 
stable against perturbations [Q. The complementary approach of the present work provides 
support to this conclusion. 



We will be concerned with two (related) quantities, the pair distribution function, g{r), 
and the static structure factor, S{k), of composite fermions. Before proceeding to the study 
of composite fermions, we discuss the behavior for non-interacting electrons, which will be 
useful when we come to composite fermions. 

The pair distribution function g{r) is equal to the probability of finding two particles at 
a distance r, normalized so that it approaches unity for large r. Two expressions for this 
quantity are 



II. NON-INTERACTING ELECTRONS 
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where the wave function \l/(ri, r2, r^v) is assumed to be normalized to unity. 
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For single-Slater-determinant (Hartree-Fock) states, g{r) can be evaluated straightfor- 
wardly by writing it as 



(3) 



where rj are single electron states and the sum is over all occupied states. For the free Fermi 
sea at zero magnetic field, this gives 



g{r) = 1 
which, for large rkp, behaves as 



kpr 



(4) 



g[r) = 1 - TT 
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For the state with n filled Landau levels, one obtains 



g(r) = 1 -e 



1 (rkpf_ 



(6) 



The Fermi sea result can be obtained from this by taking the limit n — > oo and using 

lim -Ll{-) = x-"\h{2^) . 



Fig. ^ shows how the g{r) of filled LL states evolves into that of the Fermi sea. 
The static structure factor S{k) is given by 



S{k) = -Yi{k) - N5kfl 



where 
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T{{k) =< p{k)p{-k) > 
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and <> denotes ground state expectation value. It can also be obtained from the dynamical 
structure factor S{k,u) as 



s{k) = r—s{k,uj) 

Jo n 
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S{k,uj), in turn, is related to the imaginary part of the inverse dielectric function, which, 
for non-interacting electrons, leads to 



TT 
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where n''(fc,a;) is the density-density response function. For two-dimensional spinless elec- 
trons in a Fermi sea [|1^, 
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where 



uj + ie k 
a± = — ± 
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with the square root is chosen to be on the branch with positive imaginary part. The static 



structure factor can be evaluated and found to be 
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At small fc, 



S{k) 



(16) 
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Note that there is no discontinuity at = 2fc^ in either S'(fc) or its first derivative (although 
there is one in the second derivative). 

The correlation functions gir) and S{}i) are related by Fourier transform as 



S(k) -l=p d'r e*'^"(^(r) - 1) . 



(17) 



It is also important to note that both g{r) and S{k) are properties of the ground state. 



III. COMPOSITE FERMION WAVE FUNCTIONS 



We will employ the spherical geometry in which electrons move on the two-dimensional 
surface of a sphere |T^JT^ under the influence of a radial magnetic field B originating from 
a magnetic monopole of strength Q at the center, which corresponds to a total flux of 2Q(f)o, 
where 0o = hc/e is the flux quantum. An ideal limit of zero transverse width, no disorder, 
vanishing LL mixing, and fully-polarized electrons will be assumed, allowing us to consider 
spinless electrons confined to the lowest LL moving on a disorderfree two-dimensional surface. 

The CF theory maps the interacting electrons at Q into weakly interacting fermions at 
a reduced magnetic field, corresponding to a monopole strength 

q = Q-p{N-l), (18) 

where p is an integer. The wave function of composite fermions at q (often denoted by q* 
or q'"^; the superscript will be omitted here for convenience) is related to the wave function 
of non-interacting electrons at q, denoted by $ (which is in general a linear superposition of 
Slater determinants) as 

where $i is the wave function of the lowest filled Landau level and V is the lowest LL 
projection operator. 

In general it is quite difficult to explicitly evaluate the projection operator V acting on 
the full wave function for large numbers of particles. In Ref. |^ a new and much more 
convenient approach to this problem was developed. Here, first the Jastrow factor is written 
as 

^? = I[J!, (20) 

j 

where 

r z 1 

Jj = Y[{ujVk - VjUk) exp -{(pj + 0fc) (21) 
k ^ 



and is the jth azimuthal angle. Thus can be replaced by a factor of Jj in every 
element of the jth column of the determinantal matrix. Then, the projection on lowest 
LL is accomplished by separately projecting each element. This modified wave function 
is not identical to that used previously, but very similar and has similar correlations built 
in. It should be noted that there is nothing to choose from, a priori, between different 
ways of obtaining a lowest LL wave function from the unprojected CF wave function. It is 
necessary in each case, however, to test the validity of the lowest LL wave function against 
exact solutions known from numerical diagonalization studies on finite systems. Ref. 
carried out such comparisons for systems with up to 12 particles, showing that the new 
wave functions are also extremely accurate. 

In this approach 0, ^'-^^ can be obtained from $ by replacing the single-electron eigen- 
states, Yq^n,m{^), by single-CF states, Y^j[^^{Q). Here n = 0, 1, 2, ... is the LL index, / = q+n 
is the orbital angular momentum, m = — g — n, —q — + 1, g + n labels the degenerate 
states in the nth LL, and fl represents the angular coordinates 9 and of the fermion. Y^J^^^ 
is given by 



V'^^ (O \ — AT ( -WQ+n-m + q+m q~m 

n\ / 2q + n ^ 



E(-i 
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d Y f±] J, 
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(22) 



where 



•Jj = lii^j^k - VjUk) , (23) 



k 



u = cos(^/2) exp(-i0/2) , (24) 

= sin(^/2)exp(z0/2) , (25) 

the prime denotes the condition k ^ j, and Ngnm is the normalization constant (not relevant 
for the calculation below, which involve only single-Slater determinant states). The binomial 
coefficient is to be set to zero if either (3 > a or f] < 0. 
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In order to evaluate the derivatives, we find it convenient to write them as follows: 

i d V f d V 



jp _ jp 



(26) 



where 



= Jr^J' = + > (27) 

■' OUj ^ UjVh — VjUh OUj 



Y ^ j-p ^ JP ^pY^ h— . (28) 

■' duVj ■' ^ UjVk — VjUk dvj 

For a given n, the explicit analytical form of the derivatives is used in the evaluation of the 
wave function. 

In short, Y^^^^^ can be interpreted as the single-CF wave function, from which many-CF 
wave functions are constructed in complete analogy to how many-electron wave functions are 
constructed from single electron wave functions. It should be noted that the wave function 
of a composite fermion actually depends on the coordinates of all other particles, as a result 
of the strongly correlated nature of the problem. 



IV. RESULTS FOR COMPOSITE FERMIONS 

We now compute g{r) and S{k) for composite fermions. The former is given by the 
probability of finding a pair of composite fermions at a distance r, while the latter is related 
to the equal time density-density correlation function for composite fermions. Since the 
positions and the density operator are the same for electrons and composite fermions, g{r) 
and S{k) for composite fermions are the same as those of electrons, with the expectation 
values evahiatcd in the CF wave functions. 

Two strategies will be employed for approaching the composite Fermi sea state. In one, 
we will obtain thermodynamic estimates for g(r) and S{k) for FQHE states aiu — n/ (2n-|-l) 
for n < 6, a result interesting in its own right, and look for systematic behavior leading to 
the composite Fermi sea, obtained in the limit n — > oo. The FQHE state at 
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has n filled LLs of composite fermions and occurs at effective monopole strength 



or at real monopole strength 



Q = (P+^)A^-(P+^)- (31) 

In the thermodynamic limit N oo the filling factor u = N/ (2Q) = n/{2pn + l). In various 
studies, the CF wave functions have been tested for 2/5, 3/7, 2/3, and 3/5 for systems of 
up to 12 particles, and the energies have been found to be very accurate [Elffl . We note that 



the CF state at = 1/3 is identical to the Laughlin state |]T5|, and that the non- interacting 
electron wave function $ becomes the Fermi sea state in the limit of n ^ oo. 

The second approach will be to increase the number of particles while staying at a 
vanishing effective magnetic field of composite fermions. Zero effective flux (g = 0) implies 

Q=p{N-l) (32) 

giving the compressible Fermi sea state u = l/(2p) in the limit oo. The CF wave 

functions at this flux have been studied by Jain and collaborators 0, by Rezayi and Read 
1^, and more recently by Haldane in periodic geometry [0. We will consider only the "filled 
shell" systems with N = n^, because the ground state here is a single Slater determinant, 
making computations easier, and also because it is spatially uniform (with L = 0) making 
5'(ri, a function of |ri — alone. Note that the zero-effective-fiux system with N = n'^ 
is also the smallest member of the incompressible n/{2pn + 1) FQHE state, according to 
Eq. (pOD; however, the compressible u = l/{2p) will be obtained in the limit of — oo. 

A. Pair correlation function 



Eq. (|I]) is used for a Monte Carlo evaluation of the pair-correlation function of composite 
fermions, with the distance r defined to be the arch distance between the electrons. We 
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have carried out a systematic study as a function of for FQHE states at z/ = n/{2n + 1) 
for n < 6. Fig. Q shows g{r) for the largest system studied at each fractions ||16|. There 
is no appreciable N dependence for 1/3 to 4/19, and also for 5/11 and 6/13 except for the 
last oscillation. Fig. ^ shows g{r) for g = CF states for up to 36 particles. The result 
for 9 particles is virtually indistinguishable from that obtained by exact diagonalization [Q, 
not surprising given that the CF wave function has a large overlap (0.99938 1^,^) with the 
exact ground state. 

It is illuminating to contrast the oscillations in g{r) of the FQHE and the zero-effective- 
flux states with the g{r) of a filled Fermi sea, as done in Figs. |^ and ^ An evolution into a 
Fermi-sea-like state is evident upon moving to larger values of n along the n/ {2n + 1) FQHE 
states, or upon increasing at zero effective flux, strongly indicating the presence of a Fermi- 
sea-like state at z/ = 1/2. In fact, the LLs of composite fermions bear a closer resemblance 
to a Fermi sea than do the LLs of electrons in two respects. First, more oscillations occur; 
for example, two peaks can be seen right at 1/3, the state with one filled LL of composite 
fermions, in contrast to the one filled LL of electrons, the g{r) of which does not have any 
oscillations. Second, the peak positions are much better correlated with those of the Fermi 
sea in Figs. ^ and ^ than in Fig. |I|. This implies that the formation of electron LLs is a 
stronger perturbation on the Fermi sea of electrons than the formation of CF-LLs is on the 
Fermi sea of composite fermions, or in other words, that the Fermi sea of composite fermions 
is quite robust. 

For a more quantitative analysis of the numerical results we investigate the behavior 
of g{r) at large r. For the FQHE states the oscillations are suppressed beyond a certain 
distance, as evident in Fig. |^ for 1/3, 2/5, 3/7, and 4/9, as expected from the presence of 
a gap in these systems. We are interested in the power law fall off of g{r) in the composite- 
Fermi-sea state. As seen in Fig. ^, the first few oscillations do not change substantially 
as one goes from 4/9 to 5/11 to 6/13. Therefore, we attempt to estimate the g{r) of the 
composite fermi sea by fitting the first few oscillations by the function 
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g{r) - 1 = Airkpy sin(/3rfc^ - 0) , (33) 

which provides a reasonable approximation to the actual g{r)] see Fig. |^ for a typical 
example. In all cases, we find that (3 is close to 2 (equal to 2 within the uncertainty in 
the fit), as anticipated in the above discussion. The presence of these 2kFr oscillations and 
a power-law decay of the amplitude again provides a clear evidence for the existence of a 
Fermi sea of composite fermions. The exponent a is found to be smaller than 3 but larger 
than 2 (recall that a = 3 for the non-interacting Fermi sea); its value depends somewhat on 
the range over which the fit is attempted, so should not be taken literally. 

It is natural to ask if the pair correlation function of composite fermions looks like that 
of some kind of interacting electron liquid at zero magnetic field. Interacting fermion states 
are often described by Jastrow wave functions; to make the analogy to composite fermions 
as close as possible, we try the following Jastrow wave function for interacting electrons at 
zero magnetic field: 

(^int.-FS _ JJ^ _ ^^^2(^non-int.-FS (34) 

j<k 

The pair correlation function of this state is precisely that of the unprojected composite 
fermion wave function at u = 1/2. Motivated by this, we have computed g{r) for the 
unprojected CF state at z/ = 6/13, also shown in Fig. ||. While it is qualitatively similar to 
the g{r) of the projected CF wave function, or also to the g{r) of the non-interacting Fermi 
sea, its decay is best fit by a ~ 2.4. Two conclusions may be drawn: (i) It may indeed be 
possible to find some wave function of interacting electrons at zero magnetic field, possibly 
of the form 

^int.-FS = Jl |rj - r;.|^$"°"-^'^*--^^ , (35) 
j<k 

which would behave similar to composite fermions insofar as the pair correlation function 
is concerned, (ii) The exponent a is not universal; it depends on perturbations, e.g., LL 
mixing (since the unprojected CF wave function simulates LL mixing to an extent fT^)- 
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B. Static structure factor 



In the spherical geometry, the total angular momentum L replaces the wave vector, and 
spherical harmonics are used to define Fourier transforms. The static structure factor is 
given by (for L ^ 0) 

SL = j^<\T.yL,o{Sl^)?> (36) 

3 

where the z component of L is taken to be zero without any loss of generality. The relation 
L = kR, where R = \fQl is the radius of the sphere in units of the magnetic length, will be 
used to present the results as a function of the wave vector k. The static structure function 
can be computed either directly from the wave function, or by a Fourier transform of g{r\ 
We have confirmed that both are in agreement; below we give results obtained by a Monte 
Carlo evaluation of S{}z) from the above equation. 

The static structure factors for the systems of Fig. || and ^ are plotted in Figs. |^ and 
^. Two aspects are noteworthy: S{k) are all quite similar, even for particles as few as 

= 4, and contain a sharp structure at A; = Ikp^ also reported by Haldane [0. It would 
be of interest to know if there is a real cusp at A; = 2kp. This is related to the exponent a 
characterizing the power-law fall off of g(r\ A value of a = 2 would produce a cusp in S{k\ 
as seen by noting that the Fourier transform of r"^ sin(2A;irr) is equal to 7r/2 for k < 2kF 
and arccosec(A;/2/cp') for k > 2kp. However, no cusp will occur for a > 2.0, suggesting that 
there may not be a real cusp in S{k) aX k = 2kp but only a sharp peak. The situation is 
complicated by the fact that only discrete points are available for S{k). 

The dynamic susceptibility n(A;, uo) of composite fermions has been considered within the 



framework of the Chern-Simons description in Ref. [|T^, which finds a divergent behavior as 
a; — > at A; = 2kp, but ordinary Fermi-liquid type behavior away from k = 2kp. This is 
quahtatively consistent with our results, as a divergence in S{k, oj) will produce an enhance- 
ment of S{k) at wave vectors below, but close to 2kp, while for k > 2kp, S{k) would revert 



back to Fermi-sea like behavior. However, the divergence found in |TB| is rather weak, and 
further analysis of this issue will be useful. 
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V. CONCLUSION 



Using a modified variational wave function which is easily projected onto the lowest 
Landau level, we have computed the pair correlation function and the static structure factor 
for composite fermions, and found clear evidence for Fermi-sea-like correlations in the vicinity 
of the half- filled LL, making a compelling case for the existence of a Fermi sea of composite 
fermions at = 1/2. 
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FIG. 1. Pair distribution function for non-interacting electrons at integer filling factors {v = n) 
and also for the Fermi sea. 
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FIG. 2. Pair distribution functions for several FQHE states sX v = n/{2n + 1), corresponding 
to n filled Landau levels of composite fermions. N = 50, 54, 60, 60, 60, and 54 particles have been 
used for n = 1, 2, 3, 4, 5, and 6, respectively. The statistical error in Monte Carlo is of the order 
of the noise on the curves. 
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FIG. 3. Pair distribution functions for several systems of composite fermions at zero effective 
magnetic field. The number of particles is shown on the figures in each case. 
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FIG. 6. The pair correlation function, g{r), for projected as well as the unprojected CF state 
at = 6/13. = 54 (60) particles were used for the projected (unprojected) CF wave function. 
The points are the calculated values and the curves are fits according to Eq. ( |33[ ) with a = 2.75 
(2.37), (3 = 2.11 (1.98), (p = 0.58 (-0.64), and A = 2.71 (3.63) for the projected (unprojected) CF 
state. 
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FIG. 7. Static structure factors for the systems of Fig. 0. 
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FIG. 8. Static structure factors for the systems of Fig. 0. 
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